After sequential RNase A and proteinase K treatment, DNA fragments were purified by phenol extraction and ethanol precipitation. The libraries were performed by using ThruPLEX® DNA-seq kit (Rubicon Genomics) according to manufacturer's instructions and applied to illumina platform for 151 nt pair-end sequencing by ABlife. Inc (Wuhan, China).